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o. Abstract. Overlap coincidence in a self-afiine tiling in W is equivalent to pure point 

dynamical spectrum of the tiling dynamical system. We interpret the overlap coincidence 
in the setting of substitution Delone set in R'' and find an efficient algorithm to check 
• the pure point dynamical spectrum. This algorithm is easy to implement into a computer 

I program. We give the program and apply it to several examples. In the course the proof 

of the algorithm, we show a variant of the conjecture of Urbahski (Solomyak [43]) on the 
Ifausdorff dimension of the boundaries of fractal tiles. 
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a^ ■ 

^ ■ 1. Introduction 

\ To model self-inducing structures of dynamical systems, symbolic dynamical systems 

^ ■ associated with substitutions play an important role and many works describe their spectral 

properties and geometric realizations (see |34]). To extend the symbolic substitutive systems 
to higher dimensions, self-affine tiling dynamical systems are studied in detail in |41l [26] 
- and many related studies are done along this line. These tiling dynamical systems share 

r> \ many properties with the symbolic substitutive systems and are intimately related to the 

■ explicit construction of Markov partitions. It is a subtle question to determine whether 

a given tiling dynamical system has pure point dynamical spectrum or not. It is known 
from [41 [ [26] that 'overlap coincidence' (see Def. l2.4p is an equivalent criterion to check 
this. However the overlap coincidence was not easy to compute there in practice because 
it requires topological properties of the tiles. To settle this difficulty, we shall employ the 
duality between self-affine tilings and substitutive Delone sets |22 [ I26 [ [23], An aim of this 
paper is to transfer the overlap coincidence to substitution Delone sets, find a computable 
algorithm to check the pure pointedness and implement it into a program language. 

Further motivation to show the pure pointedness comes from the study of aperiodic order. 
It is an interesting question to ask what kind of point sets, modeling atomic configurations, 
present pure point diffraction. This is related with the understanding of the fundamental 
structures of quasicrystals. It has been known from |25[ [T^ [^ that pure point diffraction 
spectrum is equivalent to pure point dynamical spectrum in quite a general setting. So the 
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algorithm we give here can be used for checking pure point diffraction of general self-affine 
quasi-periodic structures. 

There are many equivalent criteria to the pure point dynamical spectrum in literature. 
Among them, coincidences are very well known as a characterization of the pure point 
dynamical spectrum. There are many different notions of coincidences but basically they 
imply the same thing. In 1-dim substitution sequences, Dekking's coincidence is well-known 
for the case of constant-length substitutions [llj. For 1-dim irreducible Pisot substitution 
sequences or tilings, super coincidence, strong coincidence, geometric coincidence, balanced 
pairs, and boundary graph are known [TBI 121 El HOI |38l El]- In higher dimensions, modular 
coincidence was introduced for lattice substitution Delone sets [2i\ [26l I12j . and overlap 
coincidence and algebraic coincidence are known for substitution tilings and substitution 
Delone sets under the assumption of Meyer property [HI [23]. We are going to use the 
overlap coincidence for computation here. 

We note that it is essential to assume the Meyer property of the corresponding substitution 
Delone set. Otherwise, the algorithm will either not terminate as the number of overlaps 
becomes infinite, or terminate with incorrect outputs. It is shown in [27] that substitution 
Delone sets with pure point dynamical spectrum necessarily have the Meyer property. It is 
also studied in ^28j under which conditions on the expansion maps of the substitutions, the 
point sets are guaranteed to have the Meyer property. 

There are a few results in literature for the actual computation of coincidence. For 
1-dimension unit Pisot substitutions and self-affine tilings coming from their geometric re- 
alizations, computable algorithm is discussed in [371 EH] using the boundary graph. For 
irreducible 1-dimension Pisot substitution, balanced pair algorithm is implemented in [40j . 
For higher dimensions, Dekking's coincidence and modular coincidence are used for the case 
of lattice substitution Delone sets in M'^ ^[24ll26]. It was shown in [12] that the modular 
coincidence in lattice substitution Delone sets can be determined within some bounded iter- 
ations. The given bound is exponential in the number of colours of the Delone sets. It was 
conjectured in ^Jj that this bound is just polynomial in the number of colours. We give an 
affirmative answer to this conjecture not only for the lattice substitution Delone sets but 
also the substitution Delone sets with the Meyer property (see Remark 1). 

In this paper we compute overlap coincidence for general self-affine tilings. Our method 
covers, non-unit cases, higher dimensional and non-lattice based self-affine tilings. With 
regard to computation, already in the original paper by Solomyak [4jy, the number of 
overlaps becomes too large to handle by hand. Apart from 1-dimensional case with connected 
tiles (i.e. intervals), it is quite hard to check whether translated tiles have intersection. The 
implementation is already difficult for polygonal tilings, and moreover, tiles often have 
fractal boundaries in higher dimensional cases. To overcome this difficulty, we escape from 
judging interior intersection. We interpret the overlaps in terms of points and translation 
vectors, and only care distances between the corresponding translated tiles. If the distances 
are within a rough bound (see (j2.10p ). we say they are potential overlaps. Of course by this 
change, some pairs of translated tiles may not intersect, or only meet at their boundaries. 
To distinguish these cases from overlaps with interior intersection, we introduce a potential 
overlap graph with multiplicities. At the expense of having a larger graph, all computation 
becomes simple and easy to implement into computer programs. Showing that our criterion 
(Theorem 14. II (ii)) is necessary, we prove partially a variant of the conjecture which asserts 
that the boundaries of the self-affine tiles have Hausdorff dimension less than the space 
dimension d (see [l3] for the conjecture). 

The paper is organized in the following way: In Section [21 we give definitions and no- 
tations. As a main result, we present a mathematical algorithm computing the overlap 
coincidence. In Section [3] and HI we give a justification on this algorithm. In Section |5l we 
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have built a 'Mathematica' program implementing the algorithm and apply it to 1, 2 and 
3-dimensional examples. The spectral properties of some of the examples have not been 
known before. 

2. Preliminary 

The notation and terminology we use in this paper is standard. We refer the reader to |26] 
for more detailed definitions and to [21] for the standard notions. 

2.1. Tilings. We begin with a set of types (or colours) {1, . . . ,m}, which we fix once and 
for all. A tile in R'^ is defined as a pair T = (A, i) where A = supp(r) (the support of T) 
is a compact set in M"^, which is the closure of its interior, and i = l(T) S {1, . . . , m} is the 
type of T. We let g + T = [g + A, i) for g G M*^. We say that a set P of tiles is a patch if 
the number of tiles in P is finite and the tiles of P have mutually disjoint interiors. The 
support of a patch is the union of the supports of the tiles that are in it. The translate of a 
patch P hy g ^ MJ^ is g + P := {g + T : Tg P}. We say that two patches Pi and P2 are 
translationally equivalent if P2 = g+Pi for some g € W^. A tiling of M'^ is a set T of tiles such 
that M*^ = U{supp(T) : T G T} and distinct tiles have disjoint interiors. We always assume 
that any two T-tiles with the same colour are translationally equivalent (hence there are 
finitely many T-tiles up to translations). Let S(T) := {x e M"^ : 3 T,T' e T,T' = x + T}. 
We say that T has finite local complexity (FLC) if for each radius > there are only 
finitely many equivalent classes of patches whose support lies in some ball of radius R. We 
define T n A := {T G T : supp(r) n ^ / 0} for a bounded set A C W^. We say that T is 
repetitive if for every compact set K C W^, {t G M*^ : T r\ K = (t + T) fl K} is relatively 
dense. We write Bp;{y) for the closed ball of radius R centered at y and use also Br for 
Br{Q). 

2.2. Delone multi-colour sets. A multi-colour set or m- multi-colour set in is a subset 
A = Ai X • • • X Am C M"^ X • • • X M"^ (m copies) where Aj C W^. We also write A = 
(Ai,...,A„) = (Ai) 

i<m- Recall that a Delone set is a relatively dense and uniformly discrete 
subset of W^. We say that A = (Aj)j<m is a Delone multi-colour set in if each Aj is 
Delone and supp(A) := U^^Aj C is Delone. A cluster of A is, by definition, a family 
P = {Pi)i<m where Pi C Aj is finite for all i < m. The translate of a cluster P by x G M"^ 
is x + P = (x + Pi)i<m- We say that two clusters P and P' are translationally equivalent if 
P = X + P' for some x G M*^. We say that A C M'^ is a Meyer set if it is a Delone set and 
A — A is uniformly discrete (jl9j). We define FLC and repetitivity on Delone multi-colour 
sets in the same way as the corresponding properties on tilings. The types (or colours) of 
points on Delone multi-colour sets have the same concept as the colours of tiles on tilings. 

2.3. Substitutions. We say that a linear map Q : M'^ — )• M'^ is expansive if all the eigenval- 
ues of Q lie outside the closed unit disk in C. 

2.3.1. Substitutions on tilings. 

Definition 2.1. Let A = {Ti, . . . ,Tm} be a finite set of tiles in such that Tj = {Ai,i); 
we will call them prototiles. Denote by Vj^ the set of non empty patches. We say that 
Cl : A ^ V_A is a tile- substitution (or simply substitution) with an expansive map Q if there 
exist finite sets Dij C M'^ for i,j < m such that 

(2.1) n{Tj) = {n + Ti : u G Vij, i = l,...,m} 
with 

m 

(2.2) QAj = \J{Vij + Ai) for j<m. 

1=1 
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Here all sets in the right-hand side must have disjoint interiors; it is possible for some of the 
Vij to be empty. 

Note that QAj = supp{Q{Tj)) = Qsupp{Tj). The substitution (j2.ip is extended to all 
translates of prototiles by 

(2.3) n{x + Tj) = Qx + niTj), 
in particular, 

supp(il(x + Tj)) = supp{Qx + i^{Tj)) 
= Qx + Qsupp{Tj) 

(2.4) = Q(x + supp(T,)), 

and to patches and tilings by ^{P) = U{^l{T) : T G P}. The substitution Cl can be iterated, 
producing larger and larger patches Q^{P). We say that T is a substitution tiling if T is a 
tiling and f^(T) = T with some substitution J7. In this case, we also say that T is a fixed 
point of ri. We say that substitution Delone multi-colour set is primitive if the corresponding 
substitution matrix S, with Sij = ji(I'jj), is primitive. A repetitive fixed point of a primitive 
tile-substitution with FLC is called a self-affine tiling. If Q is a similarity, then the tiling 
will be called self-similar. For any self affine tiling which holds (12. 2p . we define <I> an m x m 
array for which each entry is ^ij, 

^ij = {f : X ^ Qx + d : d G Vij] 

and call $ a matrix function system (MFS) for the substitution Q,. 

2.3.2. Substitutions on Delone multi-colour sets. 
Definition 2.2. A = (A^) 

i<irn is called a substitution Delone multi-colour set in M'^ if A is a 
Delone multi-colour set and there exist an expansive map Q : and finite sets Pjj 

for i,j < m such that 

m 

(2.5) Ki = \J{QKj + Vij), i<m, 

where the unions on the right-hand side are disjoint. 

We say that a cluster P is legal if it is a translate of a subcluster of a cluster generated 
from one point of A, i.e. a + P C ^'^{x) for some k G a G M*^ and x G A. 

2.3.3. Representability of A as a tiling. Let A be a primitive substitution Delone multi- 
colour set. One can set up an adjoint system of equations 

m 

(2.6) QAj = \J{Vij + Ai), j<m 

1=1 

from the equation (j2.5p . It is known that (j2.6p always has a unique solution for which 
{^1, . . . , Am} is a family of non-empty compact sets of W^. It is proved in [22,, Th. 2.4 and 
Th. 5.5] that if A is a primitive substitution Delone multi-colour set, all the sets Ai from 
(j2.6p have non-empty interiors and, moreover, each Ai is the closure of its interior. We say 
that A is representable (by tiles) if 

A -\- A := {x -\- Ti : x G Aj, i < m} 

is a tiling of M'^, where Tj = {Ai,i), i < m, for which j4j's arise from the solution to the 
adjoint system (12. 6|) and A = {Ti : i < m}. Then A + ^ is a substitution tiling and we can 
define a tile-substitution Q, satisfying 

niA + A) = A + A 
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from ()2.6p . We call A + ^ the associated substitution tiling of A. Let <I> = {^ij) be a 
MFS for n. For any subset T = {Tj)j<m C A, ^ij{Tj) = QTj + Vij, for j < m. Let 
<I>(r) = {l->j<m^ij{^j))i<m- Then <^y(Aj) = QAj +T>ij, where z < m. For any /c € Z+ 
and X G A^ with j < m, we let <I>^(x) = ^((<I>ij(x))i<m). Note that for any /c € Z+, 
$'=(Aj) = (Q^Aj + {V%),<^ where 

ni,n2,.--i"(fe-i)<"^ 

and $'^(A) = A. 

In \22\ Lemma 3.2] it is shown that if A is a substitution Delone multi-colour set, then 
there is a finite multi-colour set (cluster) P C A for which <I>"~^(P) C <I'"(P) for n > 1 and 
A = lim„_!.oo ^"'(P)- We call such a multi-colour set P a generating set for A. 

Theorem 2.3. |26j Let A be a repetitive primitive substitution Delone multi-colour set in 
W^. Then every A-cluster is legal if and only if A is representable. 

On the other hand, if a self-affine tiling T = {Tj + Aj : j < m} is given, we get an 
associated substitution Delone multi-colour set Aj- = {A.i)i<rn of T (see [23l Lemma 5.4]). 

2.4. Pure point spectrum and overlap coincidence. Let T be a self-affine tiling in 
M*^. We define the space of tilings as the orbit closure of T under the translation action: 
Xj- = {—h + T : h ^ IR*^}, in the well-known "local topology": for a small e > two point 
sets Si, $2 are e-close if Si and 52 agree on the ball of radius around the origin, after a 
translation of size less than e. The group acts on Xj- by translations which are obviously 
homeomorphisms, and we get a topological dynamical system {Xf^W^). Let jj, be an ergodic 
invariant Borel probability measure for the dynamical system {X^,W^). We consider the 
associated group of unitary operators {Ug}g^^d on L'^{Xq-,^) : 

UgfiS) = fi-g + S). 

A vector a = (ai, . . . , a^) G is said to be an eigenvalue for the M'^-action if there exists 
an eigenfunction / G L"^ (X-j- , fJ-) , that is, / ^ and 

Ugf = e^^^f"/, for ah g G M^. 

The dynamical system (X7-, fi, M'^) is said to have pure point( or pure discrete) spectrum if the 
linear span of the eigenfunctions is dense in L'^{X-f-, fi). Recall that a topological dynamical 
system of a self-affine tiling is uniquely ergodic i.e. there is a unique invariant probability 
measure [26] . 

2.5. Overlaps. Overlap and overlap coincidence are originally defined with tiles in substi- 
tution tilings [H]. For computational reason, we define overlaps with the corresponding 
representative points of tiles here. A triple {u, y, v), with u + Ti,v + Tj G T and y G H(T), 
is called an overlap (or real overlap) if 

{u + Ai-yyn{v + Ajr^(^, 

where Ai = supp(Tj) and Aj = supp(Tj). We define {u + At — y) (1 {v + Aj) the support of 
an overlap {u,y,v) and denote it by supp(M, y, "L"). We say that two overlaps {u,y,v) and 
(«', y',v') are equivalent if there exists g such that u — y = g + u' — y' and v = g + v' , 
where u + Ti,u' + Ti G T and v + Tj,v' + Tj G T for some I < i, j < m. Denote by [(u, y, v)] 
the equivalence class of an overlap. An overlap (n, y, v) is a coincidence if 

u — y = V and u + Tj, + Tj G T for some i <m. 

Let O = {u,y,v) be an overlap in T, we define k-th inflated overlap 

^''O = {{u',Q^y,v') : u G (^^(u),v' G ^''(f), and {u\Q^y,v') is an overlap}. 
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Definition 2.4. We say that a self-affine tiling T admits an overlap coincidence if tiiere 
exists i G Z+ such that for each overlap O in T, contains a coincidence. 

Theorem 2.5. |261 [23] Let T he a self-affine tiling in such that H(T) is a Meyer set. 
Then {Xf,W^,p?j has a pure point dynamical spectrum if and only if T admits an overlap 
coincidence. 

In actual computation, it is not easy to determine whether a given triple is an overlap, 
since two points can be very close without having the interiors of the corresponding tiles 
meet. So we introduce a notion of potential overlaps. 

Let ^ G R'^ be a fixed point under the substitution such that ^ G ^{0- When there is no 
confusion, we will identify ^ with a coloured point (^, i) in Aj- with ^ G M'^. We find a basis 
of M"^ 

B = {ai, ...,ad}c H(T) 

such that 

(2.7) C + ai,...,^ + Qd G $^(0 for some £gZ+. 
Let 

Umax ■= max{|aj| : Oi G B}. 

For any n G Z-j., let 

e(") := max{|dij - d'kA : d^ G {V)ij,d'k£ G where 1 < i,j,k,£ < m}, 

Let II • II be the operator norm induced by Euclidean norm. Since Q is an expansive map, 
we can find k G such that 

IIQ"'ll<i- 

Note that for any G M"^, 

(2.8) lOM > 
Let 

(2-9) R = - 

We say that a triple {u,y,v), with u + Ti,v + Tj G T for some i,j < m and y G S(T), is 
called a potential overlap if 

(2.10) |n-y-f|<i? 
and we say that the potential overlap {u, y, v) occurs by the translation y. 
Lemma 2.6. // (n, y, v) is an overlap, then (u, y, v) is a potential overlap. 
Proof. From (j2.6p . we get 



Q^A, = l}{{V%j+A{), j<m. 

i=l 

For any i <m and a G Aj, we can write 

a = Q'^di^i + Q'^'^di^i^ +■■■ , where di„+ij„ G 
Thus for any i, j < m, a G Ai, and b £ Aj, 

\a-b\<Y^ iig-'=in<.„_, - < ^_IIq.,II - 
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If {u, y, v) is an overlap where u + Ti,v + Tj € T for some i,j < m, then 

{u + Ai-y)n{v + Aj) / 0. 
Let z G {u + Ai — y) n {v + Aj). Then z — u + y G Ai and z — v £ Aj. So 

ki — V — f < ,, ^ , ,, = R. 



Similarly to the k-th iterated overlap, for each potential overlap O = {u, y, v) in T, we 
define k-th inflated potential overlap 

$^(9 = {{u' ,Q^y,v') : u' G (^^{u),v' G <I)^(w), and {u',Q''y,v') is a potential overlap} 

and the equivalence class of ^''O 

[$^0] = {[O'] : O' G $'=0}. 

For the computation of overlap coincidence, it is important to have the Meyer property 
of S(T). The next theorem gives a criterion on Q for the Meyer property. A set of algebraic 
integers = {^i, ■ ■ ■ , Or} is a Pisot family if for any 1 < j < r, every Galois conjugate 7 of 
6j with [7 1 > 1 is contained in 0. 

Theorem 2.7. [28] Let T he a self-affine tiling in with a diagonalizable expansion map 
Q. Suppose that all the eigenvalues of Q are algebraic conjugates with the same multiplicity. 
Then H(T) is a Meyer set if and only if the set of all the eigenvalues of Q is a Pisot family. 

Summarizing the results of this paper, we provide an algorithm to determine the pure 
point spectrum of a substitution tiling dynamical system. 

Algorithm : We assume that T is a self-affine tiling in with expansion map Q for which 
H(T) is a Meyer set and Tj = i < m, are prototiles such that 

QAj = y {Vij + Ai) for j < m. 

i<m 

• Input: $ is an m X m matrix whose each entry is a set of functions from to M.'^ 
such that $ = i^ij), where = {f : x ^ Qx + d,d e V^}, i.e. $ is a MFS for T. 

• Output: True, if and only if T has pure point spectrum. 

(1) Find an initial point x such that x G ^{x). 

(2) Find a basis {ai, . . . , a^} C M.'^ such that G Uj<m((*^"'(^))i ~ for some 
n G Z+, for each 1 < k < d. 

(3) For each 1 < k < d, find all the potential overlaps Qa,o which occur from the 
translation a^- 

(4) Find all the potential overlaps Q which occur from the translations Q^a^ with 1 < 
k < d and n G lu^. 

(5) Find all the potential overlaps Qcoin which lead to coincidences within jj^-iterations. 

(6) If p{Gcoin) > piG \ Gcoin), where p{G) is the spectral radius of the graph G, return 
true. Else, return false. 
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3. Computing goincidenge 

In the rest of the paper, we assume that T is a self-affine tihng in such that H(T) is a 
Meyer set. From (j2.2p . for any Oj € Ai with i <m, we can get 

m 

(3.1) Q{Aj - aj) = [JCDij - Qaj + + {Ai - m)) for j < m. 

j=i 

We may consider new prototiles 

{Ti — ai, . . . ,Tm — Om} 

with new digit sets 

(3.2) V'ij = Vij - Qaj + ai. 

Without loss of generality we can assume that for any Tj = {Ai,i) G A, i < m, 

€ Ai. 

3.1. Meyer sets. Let A be a Meyer set and [A] be the Abelian group generated by A. 
Then [A] is finitely generated. So [A] = e-^^Zej. We define || • [| : [A] N such that 
II l^i=i Ojeill = J2i=i For each positive integer n, let 

F{n) := {u G [A] : ||n|| < n}. 

Note that F{n) is finite. Choose h > such that every open ball of radius h in M'^ meets 
at least one element in A. Since A — A := {x — y : x, y G A} is uniformly discrete from the 
Meyer property of A, we let L € be an upper bound for the number of points in A — A 
that can lie in an open ball of radius 2h. Let 

i := max{| [u| [ : u € A — A, |nl < 3h} . 

Proposition 3.1. [HI [33] Let A be a Meyer set. Then 

A-AcA + F, where F = F{2£{L-1)). 

It is proved in [26, Lemma A. 8] that the number of equivalence classes of overlaps for 
a tiling which has the Meyer property is finite. We apply the same argument to get the 
number of equivalence classes of potential overlaps for a tiling and give an explicit upper 
bound for the number. 

Lemma 3.2. Let T be a self-affine tiling and Aj- = (Aj)j<m be the associated substitution 
Delone multi-colour set ofT- Let A = Uj<m ^i- Suppose that A is a Meyer set. The number 
of equivalence classes of potential overlaps for T is less than or equal to rn^I , where 

I = #{{A + F + F + F)nBR{0)), 

with F = F{2e{L - 1)) as in Prov.[3l[ 

Proof. Let {u, y, v) be a potential overlap in T for which u+Ti,v+Tj £ T. Then \u — y — v\ < 
R. Note that u - y - v e {A - A) - {A - A). From Prop.[3TIl 

(A -A) -(A -A) c A + F-(A + F) 
C A + F + F + F . 

The equivalence classes of the potential overlaps are completely determined by i, j and the 
vector u — y — V. Thus the claim follows. □ 

Remark 1. Let A = (Aj)j<m be a substitution Delone multi-colour set for which A = 
Uj<mAj is a lattice. It has been shown in [39] that the modular coincidence, which is 
equivalent to overlap coincidence in lattice substitution Delone multi-colour sets, can be 
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determined within the exponential bound of 2'" — m — 2. Note that there are only m?! 
number of potential overlaps where I = #(A n Bfi{0)), since A — A = A. Notice that / 
is fairly small number in the case of lattice substitution Delone multi-colour set. Overlap 
coincidence of 7a can be determined within {m?I+l) number of iterations of each potential 
overlap. This polynomial iteration bound is much smaller than the exponential iteration 
bound given in [39j . 

2. We note from [27\ Th. 4.14] that A is a Meyer set if and only if H(T) is a Meyer set in 
the self-affine tiling T. 

3.2. Coincidence and computation. From now on, we assume that T is a self-affine 
tiling with an expansion map Q for which H(T) is a Meyer set. 

For a G S(T), define 

£a '■= {{u,Q"'a,v) : {u,Q'^a,v) is overlap in T, n G N }. 

For any n G Z>o, define 

Z)Qn„ :=rn(r-Q"a) 

and 

dens(Z)Q..) = ^lim —^^^^ • 

The following lemma is proved in [261 Lemma A. 9] with the subdivision graph for overlaps. 

Lemma 3.3. [2fp, Lemma A. 9] Let a G H(T). The following are equivalent: 

(i) lim„_^oo dens(I)QnQ,) = 1. 

(ii) 1 — dens{DQna) < 6r" for any n > I, for some b > and r G (0, 1). 

(iii) There exists £ G Z-|- such that for any overlap O in £a, contains a coincidence. 

The next theorem is basically in [41j and [261 Th. 4.7]. We notice here that we only need 
to consider the overlaps in £a for all a G 5 to check the overlap coincidence of T. We 
rewrite the theorem in the form that we use here. 

Theorem 3.4. [H], \2Q\ Th. 4.7] Let T he a self-affine tiling for which H(T) is a Meyer 
set. Then for any a £ 13 and any overlap O G £a, contains a coincidence for some 

(. = i{0) G N z/ and only if T admits an overlap coincidence. 

Proof. We only prove the sufficiency direction, since the other direction is clear. Suppose 
that for any a G B and any overlap O G £a, contains a coincidence. From the argument 
of [26^ Lemma A. 9], for some 6 > and r G (0, 1) 

1 — dens(L'Qn„) < 6r" for any n G N . 

Hence 

oo 

y^(l - dens(i:'Qn„)) < oo . 

n=0 

Since B forms a basis for R'^, by [41^ Thm. 6.1] the dynamical system of Delone multi- 
colour set has pure point spectrum. By Thm. 4.7 and Lemma A. 9], T admits an overlap 
coincidence. □ 

In order to find first all equivalent classes of potential overlaps which occur from the trans- 
lations of Ui for any 1 < i < d, we want to know how much region of the intersection of 7~ 
and T — cxi, 1 < i < d, we have to look. We use same notations for points with colour in 
A7- and points in M"^. This should not cause any confusion. 
Let 

e = {[P] -.P = {y,z} c Ar satisfies \y - z\ < R + \\Q-''\\amax}- 
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Let 

J{T) := {[P] E G : P C r}, where P C Ar • 
Lemma 3.5. // J(^>^(0) = Ji'^^^'^iO) for some N G Z+, then 

Proof. Let P be a cluster in <I)^+^+i(^) such that [P] G B. There must be a cluster 
p' = {y^z} C ^>^+i(^) satisfying P C $''(P'). We claim that [P'] G 9. We only need to 
show that \y — z\< R + WQ'^^Wamax- Suppose that |y — > -R + HQ"*" ||amax- Then 

11 R 

For any y' G <I>''(y) and z' G y' = Q^y + di and z' = Q^z + ^2 for some di,d2 G 

U^J<„^(P'=),^•. Thus 

\y'-z'\ = \Q^y-Q''z + di-d2\>\Q^y-Q^z\-e'^^'^ 
R 

> ii^-frii + «max - e^'') = + a„a^ from ([iJ]) 

\\Q II 

> i?+||Q~1l 

It contradicts to the choice of P. Hence [P'] G 0. From the assumption, note that 

j(f ^(0) = Ji^^'^HO) = ■■■ = Ji'^^'^HO) ■ 

So there exists a cluster P" in which is equivalent to P'. Then <1>'^(P") contains a 

cluster which is equivalent to the cluster P. Thus 

Hence 

n—>-oo 

By the repetitivity of A7-, all the clusters in A7- whose equivalent classes are in should 
occur in J(lim„_^oo Therefore = Ji^r)- □ 

Lemma 3.6. IfJ{<^>^{C)) = J(Ar), then for each aeB, 

Ga,o ■■= {[{y, a, z)] -y^ze and \y - a - z\ < R} 

contains all the different equivalence classes of potential overlaps which can occur from the 
translation of a. 

Proof. If {y, z} C A7- such that |y — a — z| < i? for a G -B, there exist n, w G Aj- such that 
y G ^^{u) and z G ^''{v). Let y = Q^u + di and z = Q^v + d2, where di,d2 G yJij<m{T^^)ij- 
Let P = {u,v}. We claim that [P] G J{<^^{C))- Suppose that \u - v\ > R + \\Q~^\\amax- 
Then 

|y-z| = \Q^u-Q^v + di-d2\>\Q''{u-v)\-e^^'^ 

1 I I (k) -Q 

— II II I ''^ ^1 6 ^ ||Q_^ II ^ CXfYiax 6 — R -\~ CXjnax 

It contradicts to the choice of {y, z} C A7-. Since J'{^'^{^)) = J{Kt), [(y, a, z)] G Ga,o- D 

Lemma 3.7. If J{^^{i)) = J{At) for some N G Z+, then for each n G Z>o and a G B, 

{[{y,Q'^a,z)] : y,z e <I)^+*''+"(^) and {y,Q"-a,z) is a potential overlap} 

contains all the different equivalence classes of potential overlaps which occur from the trans- 
lation of Q^a. 
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Proof. We argue this by induction. Note that when n = 0, the claim is true. Suppose that 
it is true for n = i, i £ Z_|_. Consider n = i + 1. Let O be a potential overlap which occurs 
from the translation of Q'^'^^a. Then there exist y,z £ A for some s G N such 

that 

= [iy,Q'+'a,z)]. 

Then there exists a potential overlap {u,a,v) with u,v such that 

{y,Q^^^a,z) G <!>{u,Q'a,v) . 

But we know that there exists a potential overlap (n', Q^a, v' ) with u',v' £ which 
is equivalent to a potential overlap {u, Q^a, v) by the assumption. Thus there exists an equiv- 
alent potential overlap {y' ,Q'^~^^a, z') to {y,Q'^~^^a, z) which is contained in ^{u' ,Q^a,v'). 
Note that 

Thus 

{[{y,Q'+'a,z)]:y,ze^''+'^+'+Hm 
contains all the different types of equivalent potential overlaps which occur from the trans- 
lation Q^'^^a. Thus the claim is proved. □ 

For any a £ B and any M G Z>o, define 

Ga,M ■= [J {[{y,Q'^0£,z)] ■.y,ze $^+'''+"(^) and {y,Q"'a,z) is a potential overlap} 

0<n<M 

Qa ■= y Ga,M and g = [J Ga- 
MeZ>o aeB 

Lemma 3.8. Let a £ B. If Qa,M = Ga,M+i for some M{= Ma) G ^>o, then 

Qa,M = Qa- 

Proof. Let [(y,Q"a,2:)] G ^«,Af+2i where < n < M -|- 2. Then there exists a potential 
overlap (y', Q^-^a, z'), with y',z' G such that 

{y,Q^a,z)£^{y',Q^-^a,z'). 

Since Ga,M = ^o,M+i; {y' ,Q"'^^ct, z') is equivalent to {y",Q^a,z") for some G 
^7V+fc+n'j^) ^j^gj.g < n' < M. So 

[<I>(y',Q"-ia,^')] = [<J>(y",Q"'«,^")]- 

Thus [(y,Q"a,z)] G^MZ+i. □ 

Let v. be the set of all equivalent classes of overlaps in T and Gcoin be the set of all 
equivalence classes of overlaps in T which lead to coincidence after some iterations. Note 
that 

Gcoin C7i C G 

and for each a £ B, Ma < m'^I, since there are at most rn^I equivalence classes of potential 
overlaps in T by the Lemma r3.21 

Lemma 3.9. // an overlap in T has a coincidence after some iterations, it should happen 
before jj^ number of iterations, i.e. 

Gco^n = {[{u,Q^a,v)] £ G : <^iiu - Q^a) n ^%{v) + % 

for some 1 < i < m and < t < \^G, where u + Ti,v + Tj £ T}. 
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Proof. Note that there are at most number of equivalence classes of potential overlaps 
in T- For any overlap O, if coincidence does not occur in $*((!?) for some < t < (jt?, 
coincidence will never occur in for any n S Z>o. Since Q is an expansive map, it is 

sufficient to check for < t < ^Q. 

4. The potential overlaps that are not real overlaps 

We aim to prove the following Theorem l4.1l in this section. The algorithm given by this 
theorem is quite simple and easy to implement and applies to all self-affine tilings whenever 
the expansion maps Q and digit sets Pjj, which define the self-affine tilings, are given. 

In the sequel, we construct a graph with multiplicities viewing potential overlaps in Q as 
vertices and define multiple edges by counting the vertices in the infiated potential overlaps. 
Hereafter we deal with the representatives of equivalence classes of potential overlaps in Q. 
Let {u, y, v) be a potential overlap, where u + Ti,v + Tj G T, Ti = {Ai,i) and Tj = {Aj,j). 
Infiating the corresponding tiles in the potential overlap (u, y, v) and intersecting them, we 
observe 

Q{u + Ai-y)r\ Q{v + Aj) 

m m 

(4.1) = [j{Ak + Vki + Qu-Qy)r\[j{A^ + V^j+Qv) 

k=i e=i 

m m 

= UU U U (.{(.^k + dki-dij + Qu-Qy-Qv)nAe) + dej + Qv). 

k=l 1=1 dki&'D^.i dijeVij 

The equivalence class [{u,y,v)] can be viewed as an element {i,y,j) where 1 < i, j < m and 
z = u — y — V with |2:| < i? where R is as defined in ()2.9p . We define the multiple edge 

(4.2) {i,z,j) ^{k,z',£) 

if z' = dki — dgj + Qxi — Qy — Qxj with \z'\ < R for dki G T^ki and dgj S T^ej, where the 
multiplicities of the edge is given by #{{dki,dij) € V^i xV^j \ z' = dki — dje + Qu — Qy — Qv} . 
Keeping the multiplicity in the graph is essential to distinguish real overlaps from potential 
overlaps that are not. Recall that {i,y,j) is a coincidence ii i = j and y = 0. We consider 
Qcoin as the induced graph of Q to the vertices which have a path leading to a coincidence. 
Also we define Gres by the induced graph generated by the complement of such set from Q, 
i.e. Gres = Q\Qcoin- For any graph G, we denote by p{G) the spectral radius of the graph 
G. 

Theorem 4.1. Let T he a self-affine tiling for which H(T) is a Meyer set. Then the 
following are equivalent; 

(i) T admits an overlap coincidence. 

(n) p{Qcoin) > p{G res I ■ 

The potential overlaps (n, y, v) can be divided into three cases. 

• No intersection overlap: (n + — y) H (t> + Aj) = 0, 

• Boundary touching overlap: u + Ai — y and v + Aj are just touching at their bound- 
aries, 

• Real overlap: {u + Ai- y)° r\{v + Aj)° / 0. 

li {u + Ai — y)r\{v + Aj) is empty, then the distance between two tiles u + Ti — y and v + Tj 
becomes larger by the iterations of the tile-substitution Vl. Therefore this potential overlap 
does not produce an infinite walk on the graph of potential overlaps by the iterations of VL. 
However, when they are touching at their boundaries, this gives infinite walks on the graph 
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and it may contribute to the number of possible paths and consequently to the spectral 
radius by repeated inflation. Our task is to prove that this contribution is small so that we 
can distinguish them from real overlaps. 

Let {V,T) be a directed graph with a vertex set V = {!,..., M} and an edge set T. 
We call {/e : e € F}, a collection of contractions /e : M"^ — )• M*^, a graph- directed iterated 
function system ( GIFS). Let F^^ be the set of edges from vertex k to i, then there are unique 
non-empty compact sets {Ek}^^i satisfying 

M 

(4.3) Ek = \J (j fe{Ee), for < M 

(see [30]). We say that (I4.3p satisfies the open set condition (OSC) if there are open sets Uk 
so that 

M 

U U feiU^CUk, for k<M 

and the left side is a disjoint union. Further if [/^ n E'fc 7^ for all k < M, then we say that 
the GIFS satisfies the strong open set condition (SOSC). 

We observe from (|4.ip 

{Ai + u - y - v) n Aj = 

(4-4) U U U U Q^' (.{(.^k + dk, - dej + Qu - Qy - Qv) n A,) + dej) . 

k<m £<m df-i&'Df^i dij&'Dij 

If {k,z',i), where z' = dki — dej + Qu — Qy — Qv, is not a potential overlap, we discard it 
from ()4.4p . 

Let M be the number of elements in Qres- Now we construct a graph for Qres identifying 
the potential overlaps in Q^es with the numbers in {1, • • • , M}. In the graph Qres, if vertices 
have no outgoing edges, we can remove them from Q^es successively. For each potential 
overlap {i,v,j) which corresponds to a vertex k < M, let 

(4.5) := {Ai + z) Aj, where z = u — y — v. 

Let be the set of edges from k to i where 1 < k,i < M. From (j4.4p . we notice that -Efc's 
satisfy GIFS (|4.3p . where fe{x) = Q~^{x + de) and de € Vgj for some £ < m. Denote by F^^ 
the set of paths of length n from k to £ and for / = ei . . . e„ G F^^ we put dj = J27=i Q^~^dei ■ 

Remark 4.2. If Qres does not contain the vertices of real overlaps, each E^ with k < M is 
a subset of dAj for some tile Tj = {Aj,j) € A. 

We use the recent development by He and Lau [15j which slightly modifies the Hausdorff 
measure. They introduced a new type of gauge function, called pseudo norm w -.W^ ^ R+ 
corresponding to Q having key properties: 

(4.6) w{Qx) = I det(g)|^/'^w;(x) 
and 

(4.7) wi^x + y) < c{w) max{'w{x),w{y)) 

for some positive constant c{w). This w induces the same topology as Euclidean norm. By 
w, they modified the definition of Hausdorff measure by: for an open set U C M'^, a subset 
X CR'^ and s,5 e M+, 

diam^(C/) = sup w{x — y), 

x,y&U 
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A C 



inf < > diam^fC/j)^ 



diam^(C/i) < J 



and 



Our new Hausdorff dimension is defined by 

dim^(X) = svip{s\ni,{X) = 00} = m{{s\ni,{X) = 0}. 
Using the pseudo norm w, one can treat self-affine attractors almost as easy as self-similar 



ones. 



To prove Theorem 14.11 we need the next lemma of Luo-Yang j29j . This generalizes a 
result in [15] and its proof basically follows from the idea of Schief [36] , but using pseudo 
norm instead of Euclidean norm. Note that strong connectedness of GIFS is essential. 

Lemma 4.3. |29|, Th. 1.1] Assume that the GIFS is strongly connected. Then following 
conditions are equivalent: 

(1) {di I I € r^^} give distinct #(r^^) points whose distance between two points has a 
uniform lower bound r > for all k,i < M and n > 1. 

(2) The GIFS satisfies strong open set condition (SOSG). 

Proof of Theorem \4-l\ (ii) =^ (i). Consider an overlap {u,y,v) where u + Ti,v + Tj E T. 
Applying to the overlapping part of the overlap, we have 

+ Ai-y)n{v + Aj))) = \ det(Q)|"^d((n + A - y) n {v + Aj)) 

where ^id is the d-dim Lebesgue measure. We know [detQI = /3 where /3 is the Perron 
Frobenius root of substitution matrix {^{T>ij)) (see ^22]). Since there are only finitely many 
overlaps up to translations, there exist r > and R > such that {u + Ai — y)° fl (u + Aj)° 
contains a ball of radius r and is surrounded by a ball of radius R. After n-iteration of 
inflation, the number of potential overlaps Kn generated from (n, y, v) is estimated: 

Cl/3" <Kn< C2/3" 

with some positive constants ci and C2. Each real overlap gives this growth of potential 
overlaps. It implies that Qres cannot contain any real overlap (Recall that we are taking 
into account the multiplicities of overlap growth). This shows the claim. 

(i) =^ (ii). We show that if all overlaps lead to a coincidence then Q^-es cannot have a 
spectral radius /3 = | deiQ\{= p{Qcoin))- By the assumption, Qres does not contain overlaps. 
So from Remark[121 Y = Ufcli ^k-, where E'fc's are defined as in ()4.5p . is the subset of the 
union of boundaries of tiles. By Lemma 14.31 the GIFS satisfies OSC because the uniform 
discreteness condition (1) of the LemmaH3] automatically follows from the fact that (X'"')ij's 
are uniformly discrete for any i,j G m and n > l(see [22j). 

We follow Mauldin- Williams [30] to compute a new Hausdorff dimension using pseudo 
norm w instead of Euclidean norm. Let 

s = dlog7/ log^, 

where 7 = p{Qres) and j3 = \deiQ\. We study the value 'H^(y). First, assuming strong 
connectedness of GIFS and OSC, we show < ?^^(y) < 00 by using standard mass distri- 
bution principle (c.f. Theorem 1.2 in [29j). Second we use a simple fact: an infinite path 
on GIFS must eventually fall into a single strongly connected component. Thus for GIFS 
without strong connectedness, we classify infinite walks on GIFS by the prefixes before they 
fall into the last strongly connected components. This gives an expression of an attractor 
of general GIFS as a countable union of contracted images of attractors which belong to 
strongly connected components. In this way we can show the Hausdorff measure is 
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positive and cj-finite, by applying Lemma 14.31 to each strongly connected component. This 
shows the new Hausdorff dimension of Y with respect to the pseudo norm w 

dim^(y) = s. 

Notice that the value ^{^{Y) can be infinite since we do not know that our GIFS is strongly 
connected (c.f. [HUl Th.4]). 

Now li s = d, then Tif^ is a translationally invariant Borel regular measure having 
positive value for any open sets because the pseudo norm is comparable with Euclidean 
norm ( [15^ Prop. 2.4]). Therefore must be a constant multiple of the d-dimensional 
Lebesgue measure, by the uniqueness of Haar measure. But this is impossible because the 
d-dimensional Lebesgue measure of the boundary of self-afiine tiles must be 0(see |32]). This 
shows s = d log 7 / log (3 < d which completes the proof. □ 

The following conjecture is a folklore. It is mentioned as an open problem in |43j from 
personal communication with M. Urbahski. 

Conjecture 4.4. For d- dimensional non-polygonal self-affine tiling T, each tile T = {A, i) 
satisfies 

d-l< dimH(«9A) < d. 

We partially solve a version of this conjecture in the following Theorem l4.5[ Indeed if the 
matrix Q gives similitudes, this settles the right inequality of Conjecture 14. 4[ 

Theorem 4.5. For d-dimensional self-affine tiling T , each tile T = {A, i) satisfies 

dim^(aA) < d. 

Proof. We consider a collection of all pairs of tiles in T whose boundaries are touching. As 
in (|4.2p and (j4.4p . we get a new GIFS which is defined on this collection. Applying the same 
argument as in Theorem l4.1l we get s = c? log 7/ log /3 < d which shows the claim. □ 

5. Examples 

We implemented Mathematica programs which perform our algorithm to check the over- 
lap coincidence for self-affine tilings. Readers can get the Mathematica programs in the 
following website. 

\protect\vrule widthOpt\protect\href {http : //mathweb . sc . niigata-u . ac . jp/\string~akiyama/Resea 

For a given expanding matrix Q and digit sets Vij of a self-affine tiling which has the 
Meyer property, the program gives outputs p{Gcoin) and p{Gres)- By Theorem 14. 1^ we can 
determine whether it satisfies overlap coincidence or not. If the tiling does not satisfy the 
Meyer property, it may not stop, or stop but produce incorrect outputs. 

In actual computation, it is the bottleneck of the program to find all initial potential 
overlaps for Lemma 13.51 and 13.61 So we use two major tricks in the program to make the 
computation fast. First, we translate the digits Vij to V[j as shown in (13. 2p such that 
the number of potential overlaps and e^'^^ are small. The size of e^^^ is significant in the 
computation of collecting all the initial potential overlaps in Lemma[331 To make e^'^^ small, 
we obtain some number of points in Ai using the tile equation (12. 2p and choose Oj among 
them which is located closest to their centroid. Then we shift tiles Ai to A^ — ai. Second, in 
order to get all the initial potential overlaps Qafi in Lemma[3]6l we try to find a fine lattice 
in R'^ such that we make an embedding of an iterated point set into the lattice taking the 
closest lattice point for each point of the set. Using the lattice, we can easily compute the 
candidates of initial potential overlaps, which is much faster than dealing with the original 
point set. We list selected examples of our computation below. 
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For 1-dimension substitution sequences, we can obtain self-similar tilings by suspension 
which associates to each letter the interval whose length is each entry of a left eigenvector 
of the incidence matrix. Pure point spectrum for the Z-action on a substitution sequence 
dynamical system is equivalent to pure point spectrum for the M-action on its suspension 
tiling dynamical system [8]. The following example shows how to obtain a tile substitution 
when a symbolic letter substitution is given. In this case, we give a separate computational 
program in the above website to check overlap coincidence directly. 

Example 5.1. We consider a period-doubling substitution — ?■ 01 and 1 — > 00. Then the 
substitution matrix is K ^ 1 and (1,1) is a left eigenvector of PF-eigenvalue 2. Giving 
length 1 to each letter, we can consider the following tile equation 

2Ai = AiU(A2 + l) 
2A2 = AiU(Ai + l) 

and get that its suspension tiling has overlap coincidence. 

We have computed a substitution 051000, 1 324100, 2 24100, 3 324333, 
4 051433, 5 51433 given in ^ Ex. 5.3] and a substitution ^ 03, 1 ^ 0, 2 ^ 21, 3 ^ 2 
[personal communication from B. Sing]. They both do not give overlap coincidence (see [39| 
Sec. 6c. 3] for alternative constructions of no pure point reducible Pisot substitution). 



We show the results of the following examples in Table[TJ 



Example 5.2. Fibonacci substitution tiling is a well-known 1-dimension tiling. The tile 
equation can be given by 



aAi 
aA2 



Ai U {A2 + 1) 



where — a — 1 
/2(x) = ax + 1. 



0. The corresponding MFS is $ 



{/i} {/i} 

{/2} 



where /i(x) = ax and 



Example 5.3. Dekking in [9l |T0] constructed self-similar tilings from endomorphisms of a 
free group. Kenyon extended this idea in [n\ §6]. We examined Example 7.5 in [41j derived 
by this method (see Fig 1(a) ). It is known that the corresponding tiling dynamical system 
is not weakly mixing and has a large discrete part in the spectrum. We check that the 
dynamical system has pure point spectrum. The tile equation is 



aA\ 
aA2 
aA-s 



A2 

{A2 

Ai- 



l-a)U{A3- 1) 



1 



with a ~ 0.341164 + il. 16154 which is a root of the polynomial + x + 1. 
We identify C with to simplify the notation; the multiplication of a in 



by the multiplication of a matrix Q 
multi-colour set is given by: 

Ai = 
A2 = 
A3 = 



a 



m 



where a = a + bi. 



is expressed 
The Delone 



aAs - 1 
aAi U (aA2 
aA2 - 1 . 



a) 
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(a) Example 15.3 



(b) Example Isn 



We take a basis B of translation vectors {I + a,a + a^} C A2 — A2 C H(T). In Tabled! we 
write this choice A2 as Colour 2. The MFS is 



{/l} {/2} 

{/s} . 

where /i = ax, f2 = ax — 1 — a and /s = ax — 1. We obtain #Qcoin = 15, #Gres = 243, 
piQcoin) ~ 1.46557 and p{Qres) ~ 1.32472. This shows overlap coincidence and therefore 
the tiling dynamical system associated with this tiling has pure point spectrum. Since this 
case is self-similar, the Hausdorff dimension w.r.t. the pseudo norm coincides with the 
usual Hausdorff dimension. So the Hausdorff dimension of the boundary of each til^ is 

2\og{p{g,o^n))/\og{p{gres)) = 1.47131. 

Example 5.4. Continuing Ex. 15. 3^ we also looked at the self-affine tiling example in [TSl 
Fig. 2 and 3] for which the eigenvalues of the expansion map satisfy — — 4x + 3 = and 
the dominant eigenvalue is not a unit. Let a ~ 2.19869, (3 ~ —1.91223 be two real roots of 



it. Set Q 



1 = (1, 1) and V = (a,/3). Then the tile equation is 



a 

fi 

QAi = (A3 + v)U(^3 + 2v-l)U(^3 + 3v-2-l) 

QA2 = {Ai + 1) U {A2 + 1) U {A2 + v) U {A2 + 2v - 1) U (^2 + 3v - 2 • 1) 
QA3 = A1VJA2 

It gives overlap coincidence as well. See Tabled] with the notation = (a^,/3^). 

Example 5.5. It is known that domino tiling does not have pure point spectrum [H 
Ex. 7.3]. The expansion map and MFS are 



Q 



2 
2 



and <1> 



{/a, 74} {/ii/e} 



i/1,/5} {/2,/4}, 

where h = Qx, f2 = Qx + (1,0), h = Qx + {0,1), /4 = + (1, 1), /s = Qx + (0,3) and 
/6 = gx + (3,0). 



"'^This number depends on other parameters we choose for computation. 

^In this example, the graph Qres is weakly connected and has only one strongly connected component of 
spectral radius greater than one. Therefore the boundary of each tile has the same dimension. The same 
holds for all examples in this paper. 
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Example 5.6. In the relation to the exphcit construction of Markov partition of toral 
automorphism, Thurston in [44j introduced {d — l)-dimensional self-similar tilings from 
numeration system based on Pisot unit of degree d which is called Pisot dual tilings. Their 
basic properties are studied in Akiyama [Ij. Such tiling dynamics are expected to be pure 
point. We confirm that Pisot dual tilings associated to x'^—x'^—x—1, x^—x—1, x^— 2x^— x— 1, 



— — 1 and x^ — 
computation. 

(1) The Pisot dual tiling associated to x^ — x — 1 has an expansive factor a 
zO. 744862 for which + a - 1 = 0. The MFS ^> is 



1 admit overlap coincidences from our algorithmic 



-0.877439 + 



V{/2} 



{/i} 







{/i} 






{/i} 





{/i} 



where /i = ax and /2 = ax + 1. 
(2) The Pisot dual tiling associated to x^ 
map Q and MFS ^ are 



■ X - 



1 is 3-dimensional. The expansion 











ai 


J_ 














J_ 








as 



and 



/{/i} 

{/2} 



V 

where 01,02,03 are roots of the polynomial x^ 
Q{x + (1, 1, 1)). The translation vectors are (2ai^- 
-a2^,2a3+a32-a3^), (2+01+01^- 



{/i} {/i} 


1/2} 



{/2} 



ai^,2a2+a2^- 



- X — 1, /i = Qx and /2 - 
ai", 2a2^ - 02^, 203^ - 03^), (2ai + 01^ - 
ai^, 2+a2+a2^-a2^, 2+03+03^-03^ 



Example 5.7. Geometric realization of 1-dimensional substitutions has been studied for 
a long time, which is motivated by Markov partition of toral automorphism. The original 
idea came from Rauzy [35] and got extended in a great deal to Pisot substitutions in [2] 
by Arnoux-Ito. They have a domain exchange structure coming from substitutions and 
inherit their spectral properties (see also [13j). Recently Arnoux-Furukado-Harriss-Ito in [3] 
generalized the idea to a special class of complex Pisot substitution which no longer has direct 
domain exchange structure but gives an explicit Markov partition of 4-toral automorphism. 
We examined the example in ^ Prop. 6.8] whose tile equation is 



oAi 
0A2 
0A3 
0A4 
0A5 
oAq 



^2 U^3 
A4UA5 

Ai 

A2 + a- 
^4 + 1- 



o 



o 



with o 
$ is 



-0.727136 + iO. 934099 which is a root of the polynomial x^ - x^ + 1. The MFS 



{/i} 
{/i} 




V 









{/l} 










{/3} 



{/l} 








{/l} 











{/l} 






^ 

* 



{/2} 



o 



and /3 



ox + o 



o 



The result shows that 



where /i = ox, /2 = ox + 1 
p{Gcoin) ^ 1.40127 and p{Qres) ~ 1.22074 and it implies overlap coincidence. We note 
that the Hausdorff dimension of the boundary of each tile is 1.18242. 
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Figure 1. Fractal Penrose tiling 



Example 5.8. Bandt-Gummelt in [5] gave two fractal Penrose tilings by fractal kites and 
darts having exact matching condition. We confirm that these tilings admit overlap coinci- 
dence. The expanding matrix Q 



matrix such that $ 



Q ^ , where _ t _ i = o and the MFS ^> is 20 x 20 
{Fij) and Fij = {f : f = Qx + d,d e Vij}. We give the digit set 

matrix (Vij) in Figuredl where w = cos(7r/5) + isin(7r/5), t = -^i^^^, and c G C satisfies 
f{t^i) = i with f{z) = jw'^ + c. Note that the Hausdorff dimension of the boundary of each 
tile is 1.26634. 

Example 5.9. Higher dimension chair tiling is discussed in [24]. We consider a 3-dim chair 









0\ 




defined by the expansion matrix 


(: 


2 


and the MFS ^> 








2/ 




/{fuh} {/l} 


{/i} {/i} 





{/i} 


{/i} {/i} \ 


{/2} {hJs} 


{/2} {/2} 


{/2} 





1/2} {/2} 


{/3} 1/3} 


{/3,/7} {/3} 


1/3} 


{/3} 


{/3} 


{M {M 


{M {/4,/8} 


{/4} 


{M 


1/4} 


{/5} 


{h} {h} {fl,h} 


{h} 


1/5} {/5} 


if 6} 


Us} {fe} 


{/s} 


{/2,/6} 


Ue} {/e} 


Ur} {fr} 


{/7} 


{M 


{M 


{/3,/7} {M 


V {fs} {/s} 


{M 


{fs} 


{/s} 


{/s} {/4,/s}/ 



where/i = Q2; + (0,0,0), /2 = gx + (1,0,0), /3 = Qx + (0, 1, 0), /4 = Qx + (1, 1, 0), /g 
Qx + (1, 1, 1), /e = Qx + (0, 1, 1), /7 = Qx + (1,0, 1), and /g = Qx + (0, 0, 1). 



Example 5.10. A 3-dimension Thue-Morse tiling can be given by the expanding matrix 
and the MFS 

Q=(0 2 0| and $ = A-^i' •^4, /e, /t} {/2, /s, /s, /s} 
Iq q 2/ \{f2, fs, fb, fs} {fi, fi, fe, fr} 

where /i, 1 < i < 8, are given as above in Ex. 15. 91 

Acknowledgments. We thank Boris Solomyak and Benoit Loridant for helpful discussions. 
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